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ABSTRACT 



Context. NGC 300 X-1 is the second extragalactic candidate, after IC 10 X-1, in the rare class of Wolf-Rayet/compact object X-ray 
binary systems exemplified in the Galaxy by Cyg X-3. From a theoretical point of view, accretion onto a black hole in a detached 
system is possible for large orbital periods only if the mass of the relativistic object is high or the velocity of the accreted wind is low. 
Aims. We analysed a 2 week SWIFT XRT light curve of NGC 300 X-1 and searched for periodicities. 

Methods. Period searches were made using Lomb-Scargle periodogram analysis. We evaluated the confidence level using Monte 
Carlo simulations. 

Results. A period of 32.8±0.4h (3cr error) was found for NGC 300 X-1 with a confidence level >99%. Furthermore, we confirm the 
high irregular variability during the high flux level, as already observed in the XMM-Newton observations of the source. A folded 
XMM-Newton light curve is shown, with a profile that is in agreement with SWIFT. The mean absorbed X-ray luminosity in the 
SWIFT observations was 1.5 x 10 38 erg s" 1 , close to the value derived from the XMM-Newton data. 

Conclusions. While Cyg X-3 has a short period of 4.8 h, the period of NGC 300 X-1 is very close to that of IC 10 X-1 (34.8±0.9h). 
These are likely orbital periods. Possibility of formation of accretion disk for such high orbital periods strongly depends on the 
terminal velocity of the Wolf-Rayet star wind and black-hole mass. While low masses are possible for wind velocities < 1000 km 
s , these increase to several tens of solar masses for velocities > 1600 km s and no accretion disk may form for terminal velocities 
larger than 1900 km s -1 . 
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1. Introduction 

Wolf-Rayet/black hole binaries are believed to be stars in the 
evolutionary stage following high-mass X-ray binaries. The ex- 
istence of helium-star/co mpact object X-ray binaries wassug- 
gested independently by Ivan den Heuvel & de Loorel d 1 973b to 
explain the nature of the g alactic source Cyg X-3, and by 
iTutukov & Yungelson! d 19731) based on results of evolutionary 
computations. To appear as an X-ray source, an accretion disk 
must form and hence the velocity of the Wolf-Rayet (WR) star 
wind must be slow enough for the mat erial around the com- 
pact o bject to be accreted. According to llllarionov & Sunvaevl 
d!975l) . a black hole appears as a strong X-ray source in a de- 
tached binary system only when the orbital period, P OI \,: 
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where Mbh is the black hole mass is solar units, viooo is the ve- 
locity of the accreted wind in units of 1000 km s~' and 5 is a 
dimensionless parameter o f or der unity. It has been shown by 
lErgma & Yungelson! d 19981) and lLommen et al.l(l2005l) that these 
periods for solar metallicity stars cannot be larger than several 
tens of hours. 

So far, Cyg X-3 is the only valid candidate in our galaxy 
for a Wolf-Rayet/compact object X-ray binary system. Its X-ray 
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luminosity is high, Lx ~ 1 38 erg s~'. The companio n star was 
identified as a WR star by Ivan Kerkwiik et all d 19921) and then 
designated as WR 145a in the 7 th catalogue of galactic Wolf- 
Rayet starsj va^e^H^h3200l]). Its orbital period is very short, 
4.8 h dParsignault et al.lll972l) . 

IC 10 X-1 (L x ~ 1.2 X 10 38 erg s -1 ), in the starburst galaxy 
IC 10 located at 0.8 Mpc, was the first extragalactic candidate fo r 
this class of objects dBauer & Brandtl 12004 IWang et al.ll2005l) . 
A period of 34.8±0.9h has been observed recently thanks to 
SWIFT observations (A. Prestwich et al. , ATel #955, paper in 
preparation). We report in this Letter the discovery of a very sim- 
ilar but slightly shorter period of 32. 8 h for NGC 300 X- 1 , which 
is the second ex tragalactic Wolf-Ray et/compact object X-ray bi- 
nary candidate dCarpano et al.1120071) . 

NGC 300 X-1 is the brightest X-ray point source in the 
dwarf spiral galaxy NGC 300 at a distance of ~ 1.88 Mpc 
dGieren et al.l 120051) . The galaxy is almost face-on and has 
a low Galactic column d ensity of Nn = 3.6 x 10 20 cm~ 2 
dDickev & Loc kman 19 93). Study of its X-ray population has 
been done by Read & Pietschl d200ll) using ROSAT and by 
ICarpano et al.l (120051) using XMM-Newton. Based on the ex- 
isting four XMM-Ne wton observations, it has been shown in 
ICarpano et al.l d2007l) that the position of the X-ray source 
(c*j2ooo = OO h 55 m 10 s :00, l W) = -37°42'1 2 , /06) coincides with a 
WR candidate, WR 41 ( Sell i Id e t al. 2003). within 071 1 + 0745. 
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WR 41 has now been spectroscopically confirmed as an early- 
type WN star (Crowther et al., in preparation). 

The four XMM-Newton light curves, lasting ~10h each, 
showed irregular variability, and during one observation, the 
flux increased by about a factor of ten in 10 h. No period be- 
tween 5 sec and 30 ksec (8.3 h) was found in the data. The mean 
observed (absorbed) luminosity in the 0.2-10keV band was 
~ 2 x 10 38 ergs~'. The unabsorbed X-ray luminosity reached 
Lo.2-iokeV ~ 1 x 10 39 erg s _1 suggesting the presence of a black 
hole, altough beamed emission from a neutron star cannot be ex- 
cluded. The spectrum could be modelled by a power-law with 
F ~ 2.45 with additional relatively weak emission, notably 
around 0.95 keV. 

In this Letter, we report the discovery of a 32.8 h period for 
NGC 300 X-l. The remainder of the Letter is organised as fol- 
lows. Section [2] briefly describes the SWIFT observations and 
data reduction. In Sect. [3] we report analysis of the SWIFT XRT 
light curve and search for periodicities using a Lomb-Scargle 
periodogram analysis. A folded XMM-Newton light curve is 
shown in Sect. |4] while a discussion of our results is given in 
Sect. |5] 




100 200 

Time (hour) 



Fig. 1. SWIFT XRT background-substracted light curve of 
NGC 300 X-l. Overplotted is the best-fit sinusoid function. 



2. Observations and data reduction 

NGC 300 X- 1 was observed w ith the SWIFT Gamma-Ray Burst 
Explorer dGehrels et al.1 12004) between 2006 December 26 and 
2007 January 10, for a total of 83 ksec. The light curve of 
NGC 300 X-l was e xtracted from the X-ray Telescope, XRT 
dBurrows et al J 120051) which operates in the 0.2-10keV energy 
band. There are 146 XRT observations lasting between 10 and 
1477 sec. We kept only data from 124 observations lasting more 
than 100 sec. 

For the production of the X-ray light curve, we analysed 
the calibrated and screened PC event files (level 2) provided 
in the set of data products. Source and background regions 
were extracted using the FT00L5Q ftselect task. The circular 
source region was centred on the XMM-Newton source position 
(anooo = 00 h 55 m 10:00, 5 J20 oo = -37°42'12'/06) with a radius 
of 40", which is larger than the telescope PSF (18"). A similar 
sized region was extracted for the background, in a blank region 
close to NGC 300 X-l. 

3. Time analysis of SWIFT light curve 

The SWIFT XRT background-subtracted light curve is shown 
in Fig. [T] Times are given in hours from the beginning of the 
observation. We overplotted the best-fit sinusoid function. For 
clarity, the amplitude has been multiplied by a factor of 1.5. It is 
clear that the flux varies in a regular way, with the minima likely 
to be eclipses of the accreting companion. 

We searched for a periodic signal bet ween 5 and 1 00 h, us- 
ing a Lomb-Scargle periodogram analysis dLombll 1 9761 : IScarglel 
Il982l) . By means of Monte Carlo simulations, we evaluated 
the confidence level assuming a null hypothesis of white noise. 
Results are plotted in Fig. [2] The full, dashed and dotted lines 
represent the 68%, 90% and 99% confidence level respectively. 
We found that the 32.84 h period is significant at a confidence 
level > 99%. To estimate the error, we fitted a sine function us- 
ing the IDL task curvefit keeping trial periods fixed. The re- 
duced chi-square, xl with v=121, is shown in Fig. [3] The full, 
dotted and dashed lines represent A% 2 =1.00, 2.71 and 6.63 re- 
spectively. The corresponding 1, 2 and 3 cr period range are 
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Fig. 2. Search for periodicities for the SWIFT XRT light curve 
of NGC 300 X-l using a Lomb-Scargle periodogram analysis. 
The full, dotted and dashed lines represent the 68%, 90% and 
99% confidence level respectively. 



[32.67-33.00], [32.57-33.12] and [32.42 -33.28] respectively. 
Note that the xl larger than 1 shows that the light curve cannot 
be described by a pure sinusoid function. 

The SWIFT XRT light curve folded at 32.84 h is shown in 
Fig. |4] Phase zero is associated to the beginning of the first 
SWIFT observation. From Fig.Q]and Fig- HI we confirm irregu- 
lar high varia bility outside the ecli pse as observed in the XMM- 
Newton data dCarpano et al.ll20071) . 

We estimate the X-ray luminosity by converting the mean 
count rate, 0.012 counts 1 to flux with WebPIMMSp , using the 
spectral parameters derived bv lCarpano et alj d2007lK The mean 
absorbed luminosity in the 0.2-10 keV energy band is 1.5 x 
10 38 erg s -1 , which is close to that fou nd from the XMM-Ne wton 
data, L .2-iokev ~ 2 x 10 38 ergs _1 dCarpano et al.1 120071) and 
close to the ROSAT value, Lo.i-2.4keV ~ 2.2 X 10 38 ergs _1 
dRead & Pietschll200ll) . 
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Fig. 3. Period error search fitting sine curve. The full, dotted and 
dashed lines represent Ay 2 =1.00, 2.71 and 6.63 respectively. 
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Fig. 4. SWIFT XRT light curve folded at 32.84 h using 30 bins. 
Phase zero is associated to the beginning of the first SWIFT ob- 
servation. 



4. XMM-Newton folded light curve 

We searched for periodicity in the four XM M-Newton data 
sample s observed between 2000 and 2005 (see ICarpano et al.l 
(I2006h for more details about the observations). Although the 
unfavourable sampling precludes a rigorous period search, we 
have used some a priori information we have on the eclipse 
profile to fold the l ight curve. From Table. 1 and Fig. 2 of 
ICarpano et al.l (120071) . we can compare the flux and light curve 
shape of the several observations to the SWIFT profile of Fig. [4] 
The flux was minimum in the first XMM-Newton observation 
and began to increase at the end. This could be associated to the 
eclipse state of the X-ray source. The second XMM-Newton ob- 
servation, 6 days later, is likely associated with an eclipse egress. 
In the third observation the flux is lower than in the second and 
fourth observations, and shows a small decrease trend: this data 
set could be associated with the beginning of the eclipse ingress. 
And in the fourth XMM-Newton observation, the flux is high 
and likely outside the eclipse. We used this information to con- 
strain the phase of each beginning of data set, in the folded light 



Fig. 5. MOS XMM-Newton light c urve folded at a per iod of 
33.066h using time bins of 300 sec (ICarpano et al.ll2007l) . First 
observation is in black, second in red, third in blue and fourth 
in green. Phase zero is associated with the start of the XMM- 
Newton observations and is by chance during eclipse as for the 
SWIFT data. 



curve. Only few values of the period around 33 h, but larger than 
32.8 h, are possible to provide a reasonable profile. 

Fig. [5] shows the XMM-Newton EPIC MOS light curve 
folded at a period of 33.066 h, within the 2 cr error of the SWIFT 
period. Phase zero is associated with the start of the XMM- 
Newton observations and is by chance during eclipse as for the 
SWIFT data. Comparing Fig. [4] and Fig. [5] it seems that both 
curves have a small dip around phase 0.6-0.8. Further observa- 
tions of the source are clearly necessary to confirm this and other 
features of the periodic light curve. 

5. Discussion 

In the evolutionary scenario for WR/compact object X-ray bina- 
ries that has been suggested to explain the short orbital period 
observed in Cyg X-3, the immediate precursor of the system is 
a neutron star or black hole orbiting an OB star. When the lat- 
ter leaves the main sequence, matter is transferred due to Roche 
lobe overflow and a common envelope forms. Due to friction, the 
distance between the early-type star core and the compact object 
decreases. A merger is avoided if the binding energy of the hy- 
drogen envelope is lower than the energy released by spiral-in, 
leading to a short-period binary system consisting of a WR star 
and a compact object. Formation of an accretion disk around a 
black hole from the strong wind of the helium star is then possi- 
ble if the orbital period satisfies Eq. (Q]i. 

We now derive possible values for the masses of the black 
hole and the Wolf-Rayet star that allow the formation of an ac- 
cretion disk for an orbital period of 32.8 h. Kepler's third law 
gives: 

(2) 



36P^ 3 (Mbh + M W r) 1/3 CR©) 



where, for a circular orbit, a is the binary separation, P or b the or- 
bital period in hours, and Mbh and Mwr the component masses 
in solar units. The velocity o f the Wolf-Rayet wind, vwr at a can 
be approximated by a yS-law dLamers & CassinelTill 1 999b : 



vwr(a) = v + (v 



*>(- "-ft 



(3) 
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Fig. 6. Allowed WR and BH masses, in solar units, that satisfy 
the condition for the formation of an accretion disk for different 
values of the terminal velocity of the wind. Different values of 
the wind velocity ft parameter (2.0, 1.5, 1.0 and 0.5) are repre- 
sented by dark to light blue. For each value of ft, the allowed 
region of Mwr-Mbh parameter space extents all the way up to 
the right. 

where Kwr is the Wolf-Rayet star radius (in R Q ), the terminal 
velocity, vo the initial velocity (~0.01voo) and the ft parameter 
describes the steepness of the law. Note that ft = 1 is the pre - 
ferred value for Wolf-Rayet winds iGrafener & Hamannll2005[) . 
Fo r the radii of Wolf-Rayet s tars, we can use the relation given 
bv lSchaerer & Maederi (1 1992b : 

log(flwR) = -0.6629 + 0.5840 log(M WR ). (4) 

Given that orbital period of the system is known, Eqs. (Q]) - (HJi 
define combinations of the component masses for which an ac- 
cretion disk may form. The plots are shown in Fig.|6]for different 
values of the ft parameter and the terminal velocity. We restrict 
to Mwr >7 M and Mbh >3 M Q , which are standard lower limit 
for the WR and BH masses respectively. 

Looking at these graphs we can note that, for a terminal ve- 
locity around 1000 km s , whatever the mass of the Wolf-Rayet 
star and the value of the ft parameter, the lower limit for the black 
hole mass is below 7 M . With higher values for the terminal ve- 
locity, this lower limit increases significantly and becomes more 



and more dependent on the ft parameter. For velocities of 1600 
km s _1 , the mass must be at least several tens of solar masses, 
while no accretion disk may form for terminal velocities signif- 
icantly higher than 1900 km s _1 . From the optical spectrum of 
WR 41, the terminal velocity of the wind is about 1250 km s~' 
and the mass of the Wolf-Rayet star is estimated between 18 M Q 
and 40 M (Crowther et al., in preparation). This leads to a black 
hole mass for NGC 300 X- 1 larger than 1 3 M for ft = 1 . 

Similar arguments apply for IC 10 X-l: the mass of the Wolf- 
Rayet star wa s derived to be 35 M a and its terminal wind velocity 
1750 km s^ 1 dClark & Crowtherl2004l) . This means that the mass 
of the black hole companion must be at least of ~35 M , for an 
orbital period of 34.8 h. 

To conclude, it seems a surprise that orbital periods found in 
both extragalactic WR/compact object X-ray binary candidates 
IC 10 X-l and NGC 300 X-l are so similar. Furthermore, their 
difference to the short period of Cyg X-3 may suggest different 
paths of evolution. 
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